The statistics of particle velocities in dense granular flows 
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We present measurements of the particle velocity distribution in the flow of granular material 
through vertical channels. Our study is confined to dense, slow flows where the material shears 
like a fluid only in thin layers adjacent to the walls, while a large core moves without continuous 
deformation, like a solid. We find the velocity distribution to be non-Gaussian, anisotropic, and 
to follow a power law at large velocities. Remarkably, the distribution is identical in the fluid- 
like and solid-like regions. The velocity variance is maximum at the core, defying predictions of 
hydrodynamic theories. We show evidence of spatially correlated motion, and propose a mechanism 
for the generation of fiuctuational motion in the absence of shear. 

PACS numbers: 45.70.Mg, 47.55.Kf 
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Granular flows are ubiquitous in nature (rock and snow 
avalanches, movement of sand dunes) and in industrial 
processes (transport of food grains, ores, pharmaceutical 
powders), yet they are poorly understood. In the regime 
of dense, slow flow there is enduring contact between par- 
ticles, and shear stresses are generated primarily by dry 
friction 0,0]- In this regime, it is often observed Q Q 
that there is fluid-like deformation only in thin layers, 
while large regions resist continuous deformation despite 
the action of a finite shear stress, like a solid. As in 
molecular fluids, knowledge of the statistics of particle 
motion will enhance our understanding of the hydrody- 
namics of granular materials, and provide insight into the 
microscopic basis of the fluid-like and solid-like response. 
Moreover, the marked difference in the way dense granu- 
lar materials and conventional fluids respond to applied 
forces prompts the question of whether there is a funda- 
mental difference in their statistical nature. 

Previous experimental studies on the statistics of gran- 
ular flows have focused mainly on the rapid flow regime, 
where grains interact through impulsive inelastic colli- 
sions [a, . Studies on the statistics of dense, slow flows 
have been far fewer, perhaps due to difficulties in probing 
dense and opaque systems. Measurement of the distribu- 
tion of the contact forces || in a monolayer of disks re- 
veal an anisotropic and inhomogencous network of "force 
chains", which has also been corroborated by particle 
dynamics simulations ||. Measurements of particle dis- 
placements by an indirect, spatially averaged technique 
[lOj and by direct particle tracking near a transparent 
wall show diffusive motion at long time scales, but 
the studies differ on their inferences of the microscopic 
mechanism at play - the former concludes that diffusion 
arises after repeated impulsive collisions between parti- 
cles, while the latter attribute it to a mechanism of rear- 
rangement in the network of abiding grain contacts. 

In our investigation, the flow of dry, cohesionless, 
spherical, almost monodisperse glass beads (Fig. lb) of 
mean diameter d p = 0.8 mm through a vertical channel 
was studied. A hopper fed the glass beads into the chan- 



nel, and the flow rate was controlled by the size of the exit 
slot at the base of the channel (Fig. la). The front and 
back faces of the channel were smooth, transparent glass 
plates, and the side walls were aluminium bars, which 
in some experiments were roughened by sticking 80 grit 
sandpaper on them. The channel width 2W and the size 
of the exit slot could be easily varied. On illuminating 
the channel from the front, the reflected light from the 
beads close to the front glass face were visible as bright 
spots (Fig. lb). Video images of the flowing beads were 
taken with a CCD camera (mounted on a translating 
stage) at 25 frames/s, for at least six minutes at each lo- 
cation in the channel, and transferred simultaneously to 
a computer. Visual inspection of the video movies clearly 
revealed long-term contact between particles, with rub- 
bing and rolling interactions predominating. Each video 
frame was analyzed, using standard image analysis tech- 
niques, to determine the ccntroids of the bright spots 
(Fig. lb) . The velocity components c x and c y of the par- 
ticles were determined from the displacement of the cen- 
troids between adjacent frames. To ensure that spots in 
adjacent frames were from the same particle, the velocity 
of a particle was determined only if its displacement was 
less than d p /2, limiting the maximum measurable veloc- 
ity to roughly lOmm/s. All measurements were made 
far enough from the entry and exit so that the flow was 
fully developed, i.e. the velocity field did not vary in the 
direction of flow. 

The profiles of mean properties were determined by no- 
tionally dividing the channel into vertical bins of width 
dp and averaging the relevant particle properties in each 
bin over the entire duration of the video movie. Figure lc 
shows a typical result for the mean velocity u x and the 
root-mean-square (rms) velocity fluctuation v. The mean 
velocity increases sharply within the shear layer, from a 
small but finite slip at the wall to the maximum value uq, 
and remains constant at uq in the substantial core. This 
feature of shear banding is in agreement with previous 
work on flow in channels 0, and cylindrical Couette 
cells 0, . The dependence of the shear layer thickness 
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FIG. 1: (color online) (a) A sketch (not to scale) showing the 
construction and dimensions of the channel, (b) A snapshot of 
the flowing glass beads, showing the bright spots of reflected 
light, (c) Profiles of the mean velocity u x , and the root-mean- 
square velocity fluctuation v for a rough walled channel with 
W7d p = 31.7. 



on W and the material and wall properties, and a com- 
parison of the observed mean velocity field with predic- 
tions of available hydrodynamic theories will be reported 
elsewhere 01 . Though the shear rate is negligibly small 
in the core, we find that v is highest there (Fig. lc); this 
is a surprising observation, as we expect the conversion of 
energy from mean to fluctuational motion to be high in 
regions where the shear rate is high. When a liquid (or a 
granular medium in the rapid flow regime |l4j l is sheared, 
the shear rate is high in regions where the temperature 
is high, as its viscosity decreases with rise in the temper- 
ature. Our observation of the opposite here illustrates 
the difference in the micromechanics between dense fric- 
tional flows and molecular or rapid granular flows: the 
velocity fluctuations here arise from complex, collective 
motion, with particles contacts being sustained almost 
continuously, and not from particles flying between col- 
lisions with high velocity. 

The distribution of velocities in the x and y direc- 
tions (vertical and horizontal directions, respectively, see 
Fig. la) were determined as a function of position in the 
channel for a range of the channel width, with rough and 
smooth side walls. Our visualization technique could not 
resolve movement in the z direction. The probability 
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FIG. 2: (color online) The probability distribution function 
for vertical (a) and horizontal (b) particle velocity in the fluid- 
like shear layer. The data are for four different channel widths 
with rough walls (see text), the legend giving the ratio W/dp. 
The unfilled symbols in the main panels represent data for 
negative velocity fluctuations, reflected onto the positive side 
for accommodation in the log-log plot; the filled symbols rep- 
resent data for positive velocity fluctuations. The inset in 
each panel shows the distribution in linear-log axes. The dot- 
ted line in the inset is the Gaussian distribution with the same 
variance. The straight lines in the main panels are power-law 
fits for large velocity magnitude. 



distribution function f x (£x), where £ x = (c x — u x )/v is 
the dimcnsionless velocity fluctuation in the vertical di- 
rection, was determined by making a histogram of the 
number distribution of £ x , and normalizing it so that 
ffx(£x)d4x = 1; the distribution function f y (£ y ) for hor- 
izontal velocities was determined similarly. Figure 2 
shows the distribution function in the shear layer for 
rough walled channels of four different channel widths. 
We note that f x is not symmetric about £ x = (Fig. 2a), 
as gravity induces a preference for downward velocities. 
Its decay is more rapid for negative £ x than for positive. 
There is no preferred direction for horizontal velocities 
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and the distribution is therefore symmetric about £ y = 
(Fig. 2b). In both directions, the distribution function 
decays as a power law, /,•(&) ~ when is suffi- 

ciently large. 

A striking feature of Fig. 2 is the collapse of the data for 
all the channel widths into a single curve. We find that 
data at other locations in the shear layer and for smooth 
walls (not shown) also collapse onto the same curve. As 
the flow rate, and therefore the shear rate, varied between 
all these experiments, this implies that f x and f y do not 
depend on the shear rate, but only on the local value of 
v. This result is further reinforced when we compare the 
data in the shear layer and the core (Fig. 3); the distri- 
bution function is identical in the two regions. This is 
a remarkable result, as the nature of the flow in the two 
regions is very different - there is continuous deformation 
in the shear layer, but none in the core. Curiously, this 
is similar to molecular fluids at thermodynamic equilib- 
rium, where the Gaussian distribution holds regardless 
of the state of aggregation. However, the distribution 
function away from equilibrium in molecular fluids de- 
pends on gradients of the velocity and temperature 
Granular flows are fundamentally far from equilibrium 
(indeed the flowing state bears no resemblance to the 
equilibrium, static, state), and hence the shear-rate in- 
dependence is a significant deviation from the statistical 
behavior of classical fluids. 

Our results are consistent with those of Ref . [JjJ , who 
observed anisotropy and non-Gaussian distribution of the 
particle displacement over a time scale of 1 ms. How- 
ever, they found the distribution function changing to a 
Gaussian on coarse graining in time, but we find it dis- 
tinctly non-Gaussian even with though our camera speed 
was 40 times slower. We found no perceptible change in 
the distribution function when a faster video camera (60 
frames/s) was used 0, suggesting that our data is ro- 
bust and valid over a range of time scales. 

The power-law exponents were determined by fitting 
a straight line through the data in the log-log plots for 
large velocities (Figs 2 and 3). For the horizontal ve- 
locity the exponent n y is 3.6, and for the vertical ve- 
locity the exponents are n x =2.9 for positive £ x and 
n~ = 7.4 for negative £ x . The faster decay of the dis- 
tribution for negative £ x is a reflection of the preference 
for positive (downward) velocity fluctuations that grav- 
ity induces. We do not, however, expect this form of 
the distribution to persist to very large velocities for the 
following reason: if we consider the mean velocity fluctu- 
ation square, or the specific fluctuational kinetic energy 
E3 V L = JZl c * - u *) 2 f(°x)dc x + f^c* f{cy)dc y , we 
find that the first integral is unbounded because n+ < 3. 
Thus, for the mean kinetic energy of particles to be fi- 
nite, there has to be a faster decay of f x for large positive 
c x . This can be verified if higher particle velocities are 
measured, perhaps by high speed video imaging [fillllj. 

The presence of velocity fluctuations even in the ab- 
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FIG. 3: (color online) Comparison of the probability distri- 
bution function in the shear layer with that in the core. The 
layout and symbols are as in Fig. 2. The suffix 's' and 'c' 
in the legend indicate data for the shear layer and the core, 
respectively. 



sence of a shear rate has been established in earlier stud- 
ies 0, on dense granular flows. Our observation 
that v is highest in the core is somewhat more sur- 
prising, as hydrodynamic theories 0] give the trans- 
fer of kinetic energy from mean to fluctuational motion 
as the product of the shear stress and the shear rate, 
predicting that v is minimum at the core. A plausible 
explanation for the production of velocity fluctuations 
in the absence of shear is provided by considering the 
spatial correlation of the velocity fluctuation, defined as 
Cy{Y) = (£(t,y)Z(t,y + Y)). Here ~%(t,y) is the in- 
stantaneous fluctuation velocity averaged over particles 
in the bin centered at the reference point y, c' x (t, y + Y) 
is the same in the bin centered at y + Y , and the angle 
brackets indicate an average over time. Figure 4a shows 
that C y (Y) decays rapidly with Y, in the shear layer 
and the core, but not to zero. There is a finite corre- 
lation at large distances, implying a correlated solid-like 
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FIG. 4: (color online) Correlation functions of the veloc- 
ity fluctuation in the shear layer and the core (see text for 
definitions of C y and Ct) for a rough walled channel with 
W/dp — 31.7. The reference point y was 4d p from the left 
wall for the shear layer, and channel center for the core, (a) 
The spatial correlation as a function of the horizontal separa- 
tion Y, and (b) the time correlation as a function of the time 
delay T. 

motion of the particles superimposed over the uncorre- 
cted fluid-like motion. We believe that this solid-like 
motion is a result of the stick-slip 0,0], or jamming- 
unjamming, phenomena in granular materials. When the 
medium is released, or unjammcd, during the slip phase, 
the unbalanced force accelerates the particles for a short 
period of time, and a part of the kinetic energy acquired 
by the particles is retained in fluctuational modes. The 
energy is lost in the next stick or jammed phase, but the 
cycle is repeated regularly. The decay of the time corre- 
lation, C t {T) = (4(t,y)4(i + T,y)), to zero in Fig. 4b 
indicates that the stick-slip motion itself is uncorrelated 
beyond a time scale of roughly half a second. (Note that 
Ct (T) is not the usual autocorrelation function, as we are 
not tracking the velocity of individual particles in time). 
The relative unimportance of shear in generating fluctu- 
ational motion implies that shear work largely dissipates 
as true heat through rubbing and rolling friction 1 1 7I | . It 
appears likely that the jamming-unjamming cycle is af- 
fected by shear, but it is not clear why shear diminishes 
v but leaves the velocity distribution unchanged. 

The above discussion leads to the important question 
of whether the velocity distribution we have reported is 
universal for slow granular flows, or is limited to flows 
driven by a body force. This is an issue of considerable 



importance, and can only be answered by conducting ex- 
periments in other geometries. A kinematic variable that 
we have not measured in this study is the rotation veloc- 
ity of the particles. Recent studies that have modelled 
granular media as Cosserat continua [2(3, El have shown 
the rotation velocity to play an important role in deter- 
mining the kinematics. Though its measurement is more 
difficult, it will be useful in critically evaluating hydro- 
dynamic theories and in understanding the microscopic 
difference between the solid-like and fluid- like regions. 
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